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We theoretically study the propagation of light in a Fermi- 
Dirac gas in the presence of a superfluid state. BCS pairing 
between atoms in different hyperfine levels may significantly 
increase the optical linewidth and line shift of a quantum de- 
generate Fermi-Dirac gas and introduce a local-field correc- 
tion that, under certain conditions, dramatically dominates 
over the Lorentz-Lorenz shift. These optical properties could 
possibly unambiguously sign the presence of the superfluid 
state and determine the value of the BCS order parameter. 
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After the first observations of atomic Bose-Einstein 
condensates |l| there has been an increasing interest in 
the studies of Fermi-Dirac (FD) gases [||-||- One espe- 
cially fascinating property of FD gases is that with ef- 
fectively attractive interaction between different particles 
the ground state of the system may become unstable with 
respect to formation of bound pairs of quasi-particles or 
the Cooper pairs ||lO| , [ll| . This effect is analogous to the 
BCS transition in superconductors. The particles near 
the Fermi surface having opposite momenta and differ- 
ent internal quantum number tend to appear in pairs. 
This leads, e.g., to a finite energy gap in the excitation 
spectrum of the system and to a nonvanishing anomalous 
expectation value (■!/'|(r)'0x(r)). Here the two internal 
sublevels are referred to as t and J,. 

In this paper we study the optical response of a su- 
perfiuid state in a zero-temperature FD gas at low atom 
densities for low- intensity light. We show that the pres- 
ence of the coherent pairing between atoms in different 
hyperfine levels may dramatically enhance optical inter- 
actions and the scattering of light in FD gas. 

One particularly promising candidate to undergo the 
BCS transition and to become a superfluid is spin- 
polarized atomic ^Li. Atoms in two different internal 
levels can interact via s-wave scattering and the ^Li atom 
has an anomalously large and negative s-wave scattering 
length a 2± — 2160ao, where uq is the Bohr radius. The 
nuclear spin states nii = 1 and of ^Li have been pre- 
dicted to undergo a superfluid transition at 10~^ K with 
the density of lO^^ cm^^ |,|. 

We study the propagation of light by introducing the 
dipole approximation for atoms and the corresponding 
Hamiltonian in the length gauge obtained in the Power- 
Zienau-WooUey transformation ||l^. FD gas is assumed 
to occupy two different hyperfine levels |g, T) and \g, i) 
of the same atom. For simplicity, we consider here a 
situation where there are only two electronically excited 



levels |e,t) and |e, j). All the dipole matrix elements 
for optical transitions between the levels are assumed to 
vanish, except dg|eT and dgi^i, where dgfet denotes the 
dipole matrix element for the transition |e,t} — > |5,T}- 
For instance, the electric dipole moment between levels 
t and J, vanish, if the levels refer to different nuclear 
spin states that are decoupled from the levels involved in 
optical transitions by an external magnetic field. 

In the absence of the driving light field atoms in the 
electronic ground state are described in second quantiza- 
tion by the Hamiltonian density Tii |l(]| , [ll|] : 

(1) 

where ipg^{r) is the atom field operator for level \g, v) in 
the Heisenberg picture, ^g^ is the corresponding chemi- 
cal potential, and H^^^ denotes the center-of-mass (cm.) 
Hamiltonian. We have approximated the finite-range in- 
terparticle potential by a contact interaction with the 
strenght given by Ug — Airagh/m, where ag is the s- 
wave scattering length and m is the mass of the atom. 
The atoms in different hyperfine levels can interact via 
s-wave scattering. On the other hand, there only is a 
very weak p-wave scattering between two atoms in the 
same hyperfine state which is ignored in Eq. ([l|). 

The driving light field introduces additional terms for 
the system Hamiltonian. In the length gauge the ba- 
sic dynamical degree of freedom for the light field is the 
electric displacement D(r) that interacts with the atomic 
polarization P(r) 

H2 = --P(r).D(r). (2) 

eo 

In the present case the positive frequency component of 
the polarization is given by 



(3) 



"(r) = dg|eTV'g|(r)V'eT(r) +dgXeiV'gx(i-)V'ei(r) 
^P+(r)+P+(r). 

The polarization self-energy was shown in Ref. ||l^ to 
be inconsequential for dipole atoms. We assume that 
to leading order all remaining interactions between the 
atoms and the light field, that cannot be accounted for 
when the atoms arc modeled as point dipoles, are gov- 
erned by the following interactions: 

u 

+ X! ^'^g'^ea^Ull^L^eal/jgi^ ■ (4) 



Here Ug^ea = '^T^'hagvea/'rn and Ue = 'iTrhae/m describe 
the two-body s-wave scattering between the atoms. For 
simphcity, the frequency of the optical transition uq is 
assumed to be independent of the hyperfine level. For 
typical values of the optical linewidth the cm. motion 
for the excited atoms may be omitted 0. 

The positive frequency component of the electric field 
E"*" may be obtained by solving the Heisenberg equations 
of motion Q 

eoE+(r) = D+(r) + J d\' G(r - r') P+(r') , (5a) 
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Here D^ is the positive frequency component of the 
driving electric displacement with the frequency f2, and 
k = Q/c. The monochromatic dipole radiation kernel 
G(r) is precisely the classical expression of the dipolar 
field, including the delta function at the origin p^ . 

In the limit of low light intensity we have derived from 
the Heisenberg equations of motion a hierarchy of equa- 
tions for correlation functions involving atomic polariza- 
tion and atom density |lj,|l^ . In the case of the present 
system we may proceed similarly. As far as the optical 
response is concerned it is again assumed that we can 
concentrate on the dynamics of internal degrees of free- 
dom for the atoms and the light. Hence, in the equation 
of motion for the atomic polarization the kinetic energy 
of the atoms is neglected. 

Light has to be present in order to produce population 
in the electronically excited levels. In the present paper 
we consider the limit of low intensity of the driving light. 
This is done by retaining only those products of operators 
that involve at most one excited state field operator or 
the driving electric displacement [H. Then, e.g., the 
term proportional to Ue in Eq. (^ has no contribution to 
the equation of motion for P+(r). 



To simplify further we assume that d 
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is useful to introduce the following projection operator 



IgTeT^lgTeT/IClgTeTl ' 



(6) 



that projects the internal degrees of freedom onto the 
subspace defined by the four hyperfine levels in consid- 
eration. For the expectation values we use the notation 
Ply = (P+), for V denoting the hyperfine state. The 
steady-state solution of Pii, is given by 

Pi,.(ri) == ap^P ■ 'Dj,{ri) + ^ J^CTyP2(ricr; rii^) 

+c^Yj I '^^'■2 P • G'(ri - r2)P2(rii^; raa) . (7) 

Here a — —V^ /[Tieo{S+ij)] is the polarizability of a single 
atom, V is the reduced dipole matrix element, p^ denotes 
the atom density in level i/, j = V^ k^ / {Gnheo) is the 



spontaneous linewidth, and S the atom-light detuning. 
We have also defined 

P2(rii.; r^a) s (^t^(ri)P+(r2)^,.(ri)), (8) 
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The normally ordered expectation value P2(riz^; r2(T) de- 
scribes correlations between an atomic dipole at r2 in 
hyperfine level a and a ground state atom at ri in hyper- 
fine level v. The tensor .7^^ generates the coUisionally- 
induced level shifts. 

Due to the hard-core interatomic potential we remove 
the contact dipole-dipole interactions between different 
atoms. In Eq. (|7|) this is done by introducing the prop- 
agator G^j(r) = G,;j(r) -f 5ijS{r)/3. The purpose of this 
definition is to yield a vanishing integral for G'(r) over 
an infinitesimal volume enclosing the origin pSl . 

So far, we have obtained a steady-state solution for 
the atomic polarization (|^) that acts as a source for the 
secondary radiation in Eq. (pa). Equation (M) involves 
unknown correlation function P2. Basically, we could 
continue the derivation and obtain the equations of mo- 
tion for P2 and for the higher order correlation functions. 
This would eventually result in an infinite hierarchy of 
equations analogous to the equations in Ref. [Q. How- 
ever, even in the case of a simple level structure and in 
the absence of the s-wave interactions the solution for the 
whole system by stochastic simulations is demanding on 
computer time H . In the studies of the refractive index 
of a quantum degenerate Bose-Einstein gas Morice et al. 
p6[ derived a density expansion in terms of the number 
of atoms repeatedly exchanging a photon by introducing 
certain approximations in the ground state atom corre- 
lations. Although the lowest order density correction for 
the electric susceptibility of a zero-temperature ED gas 
may be obtained analytically ||9[] , in the presence of non- 
trivial statistical position correlations a rigorous density 
expansion is in most cases a very challenging task. In 
this paper we consider low atom densities (in terms of 
p/k^) and approximate Eq. (0) by the decoupling that 
corresponds to the lowest order correction in Ref. ||I^ 

P2(rii/;r2cr) ~ p2{rii^,r2a)Pia{r2)/pa , (10) 

where the ground state pair correlation function p2 is 
defined by 

P2(rii^,r2cr) = (i/;t^(ri)V'J^(r2)V'g^(r2)Vv(i"i)) • (H) 

It is important to point out that the predictions of the 
expansion by Morice et al. [^ were tested for a zero- 
temperature FD gas in one dimension [g|. The agree- 
ment with the exact solution obtained by the numerical 
simulations was found to be semi-quantitative and in the 
low-density limit excellent. 

Before the light is switched on, the system is described 
by the Hamiltonian density H = Hi [Eq. (|l|)]. The as- 
sumption that the driving light only weakly disturbs the 



system allows us to evaluate the pair correlation functions 
for the ground state atoms [Eq. (|ll|)] from Hi even in the 
presence of the driving light. We assume a homogeneous 
sample and introduce a plane wave basis for the field op- 
erators: 4'gu{^) = ^~^^^X]k^ki/exp(ik-r). The Hamilto- 
nian (|l|) is diagonalized by the standard canonical trans- 
formation to the Bogoliubov quasi-particles pO| , [ll[ 



X- 



ttk = Uk^ki -i^k^-kT' 
/3-k = Wk^-kT + f'k^k I , 



(12a) 
(12b) 



where u^ and wj. are real, depend only only on |k|, and 
satisfy u^ -t- u^ = 1. The requirement that Tii in Eq. (|l|) 
is diagonal in the quasi-particle representation sets an 
additional constraint and we obtain 



ul = l\l 



Sk 



1 



a 



--2l^-f 1' (^^) 



where E^, = ^AM^, Ck = e^,- ^- + TiUg{p^+ pi)/2, and 
the energy gap A = ~fiUgV~^ I]k'"k«k(l - n-ak - S/3k)- 
In equilibrium, the quasi-particle occupation numbers 
rick = (o^kO^k) and np\^ = {Pl_(3-k) satisfy FD statistics 
with flak — npk — (e^''/'^^^ -I- 1)^^. The dispersion re- 
lation for free particles is given by ek — h^k^ /{2m) and 
the average of the chemical potentials is /I = (/i| +/i|)/2. 
For simplicity, from now on we assume /i| = l^-i- 

In the superfluid phase transition the atoms in the dif- 
ferent hyperfine levels t and I form a quasi-particle con- 
densate that results in a nonvanishing anomalous cor- 
relation {ipi{ri)ipi{r2))- The effect of this macroscopic 
two-particle coherence on the pair correlation function 
( pT| ) can be clearly seen by considering the ground state 
of Til [Eq. (0)] that is the vacuum of the Bogoliubov 
quasi-particles [Eq. (12)]. Then (for v ^ a) 



P2(vxV,V2(j) = PuPa + |(V'gi.(ri)'0gcr(r2))|^, (14a) 

P2{viV,V2v) ^ pl~\{i;l{v^)^g,{v2))\' . (14b) 

The optical response may now be evaluated by elim- 
inating D;^ and P2 from Eqs. (||), (0), and (^. Be- 
cause we are dealing with a linear theory, the electric field 
and the polarization are related by the susceptibility as 
P+ = eoxE+. We consider a situation where FD gas fills 
the half-infite space z > Q. We assume that the constant 
atom densities for the hyperfine states are equal p^ = 
Pl = p. For simplicity, it is also assumed that the s-wave 
scattering length Ogi/eo- is independent of hyperfine states 
u and a resulting in :F|| = J^n and T^i — J^ip The in- 
coming free field is written Di?(r) — Dp eexp (ikz), and 
we assume that the polarization e is parallel to the elec- 
tric dipole moments dg|eT and dgi^i- With the ansatz 
Pit(r) = Pu(r) = PeexTp{ik'z) for Im(fc') > 0, by us- 
ing Eq. (|lj), and by ignoring the effects of the surface of 
the atomic gas B] , we obtain a spatially constant suscep- 
tibility for the sample as 



A;2 



2ap 



1 - 2ap/3 -I- El + S2 



(15) 



with 



--/dVe-^-G'(r)[|(V,T(r)^.i(0))P 



S2 = 



-|(V4(r)^.i(0))| = 
--^Ti,aP2{ri^,ra). 



(16) 
(17) 



Here we have used the obvious relations /92(ri 1, V2 [) = 
P2(ri i,r2 T) and p2(ri T,r2 T) =P2(ri i,r2 [). 

In an uncorrelated atomic sample the atomic posi- 
tions are statistically independent and the pair correla- 
tion function ( pT|) satisfies p2{vv,v'a) — pvPa- In this 
case, and in the absence of the s-wave scattering, we 
would obtain Eq. ( p^ with Si = E2 = 0. This is the 
standard column density result stating that susceptibil- 
ity equals polarizability of an atom times atom density. 
Equation (|l^) also involves the Lorcntz-Lorenz local-field 
correction in the denominator. 

The quantum statistical corrections to the column den- 
sity result are introduced by Si. It describes the modi- 
fications of the optical interactions between neighboring 
atoms due to the position correlations. The second term 
in Eq. (16) represents the quantum statistical contribu- 
tion to the scattering process in which a photon emitted 
by an atom in hyperfine level v at position r is reabsorbed 
by another atom in hyperfine level v and located at the 
origin. According to FD statistics two fermions with the 
same quantum numbers repel each other and FD statis- 
tics forces a regular spacing between the atoms. The 
optical interactions are dominantly generated at small 
interatomic distances and the corrections to the suscep- 
tibility due to the second term in Eq. dlq) correspond 
to inhibited light scattering. In the absence of superfluid 
state FD gas exhibits a dramatic line narrowing Wm ■ 

The first term in Eq. (O) represents the quantum sta- 
tistical corrections to the reabsorption process between 
atoms in different hyperfine levels due to the two-particle 
coherence. This term is nonvanishing only in the pres- 
ence of a superfluid state. Because the atom pairs in 
Eq. M) interact in the triplet s-state and the total spin 
of the pair is an integer, the pairs behave as bosons [ pX)[ . 
According to the Bose-Einstein statistics two bosons at- 
tract each other and the presence of the BCS pairing 
favors small interatomic spacing, hence, enhancing the 
optical interactions and the light scattering. 

The line shift induced by S2 [Eq. lK%] is generated by 
the s-wave interactions. As far as they can be considered 
local on the scale of the optical wavelength in Eq. (Q) 
the collisions induce a local-field shift analogous to the 
Lorentz-Lorenz shift. 

The optical line shift for the atomic sample is obtained 
from Eq. (fsl) 



5/7 = 47rp + {ug - Uge)p2(T, i)/p - 67rRe 



(18) 



where we have dropped the equal position coordinates in 
P2, used p2{'ri>,'ri') — 0, and defined the dimensionless 
variables p — p/k^, P2 = P2/k^, S = S/j, a = —6n/{6 + 
i), and u^ = U(^k^ /j. The first two terms form the local- 
field shift. For ^Li the local-field shift due to the s-wave 
scattering in Eq. ( Jig ) is larger than the Lorentz-Lorenz 
shift, if 7 < UO[l + {A/nugp)'^]{ag-age)/{aoX^) pm^s-\ 
where A is the optical wavelength. Because {A/hugp)'^ is 
expected to be of the order of one Q , the local-field shift 
could strikingly depend on the BCS order parameter A. 

If the the effective range r,j of the triplet s-wave po- 
tential in Eq. (0) is very short r„ ^ 1/fc, the resonant 
dipole-dipole interactions may suppress the effect of the 
s-wave scattering on the line shift just as they cancel the 
effect of the polarization self-energy m^. However, for a 
metastable state 7"^ may be large on the time-scale of 
the atomic interactions. In that case the coUisional shift 
could be observable even for very small r„. 

To calculate the nonlocal linewidth and line shift from 



integral (16) we need to evaluate the spatial correlation 
functions by using Eqs. (Ill) and (p^). For instance, the 
anomalous expectation value reads 



(V'i(r)V'T(O)) 









riak - "-/3k) • (19) 



radiation already involves a high-frequency cut-off |1J] 
that regularizes small r behaviour. We have also plot- 
ted the case with the cut-off fee = l/?'« (the dotted line) 
with the realistic value r„ = lOOao of the triplet s-wave 
potential B. We found that the linewidth is almost in- 
dependent of the cut-off from r^ = SOOag to r„ = 0. The 
line shift from the unregularized anomalous correlation 
in integral (|l^) diverges logarithmically for small r. Al- 
though the radiation kernel (5b) involves a cut-off, the 
Lamb shift is not treated rigorously. However, for the 
present purposes we may at least obtain an estimate for 
the shift by using the T-matrix or the cut-off kc — l/r„ 
as in the case of the linewidth. For the BCS state, even 



with u„ = 



'ffe 



and moderately small 



shift is increased. 
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The chemical potential is solved from p^ ~ Pvifi)- Here 
('0x(O)V't(O)) = ~^l{'^'^g) is ultraviolet-divergent result- 
ing from the assumption of the contact two-body inter- 
action in Eq. (H^). This interaction is momentum inde- 
pendent and is not valid at high energies. To estimate 
the pairing we remove the high-energy divergence by re- 
placing hug by the two-body T-matrix obtained from the 
Lippmann-Schwinger equation y. This is done by sub- 
tracting from Eq. ( |l9| ) ^ exp(ik-r)A/(2^kl^)- As argued 
in Ref. [^ the use of the T-matrix may seriously under- 
estimate the overlap in the case of large scattering length 
\o-g\ 3> Ttj. Hence, we calculate Eq. ( |l9| ) also by introduc- 
ing a high- momentum cut-off kc- 



We have plotted the line shift from Eq. (18) and the 
linewidth r/7 = 1 — 67rIni(Ei/Q:) by assuming, for sim- 



plicity, Ug = 



J-ge, 



A — 900 nm, and a moderate value 



= — 1200ao. For the gap parameter at T = we use 



the weak coupling approximation A ~ l.TGfcsTc |l( 
where 



kuTr. 



8eF ^-2 
e ' exp 

■K 



2kF\a\ 



(20) 



with 7 ~ 0.5772 and kp = (67r2p)i/3, 

In Fig. n (a) the solid line represents the linewidth in 
the absence of the superffuid state (A = 0). The line 
narrows as a function of the density H . The presence of 
the superfiuid state broadens the line. The linewidth is 
finite even without the regularization in the anomalous 
correlation (the dashed fine). This is because the dipole 
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FIG. 1. The optical (a) linewidth and (b) line shift as a 
function of the atom density per cubic optical wave number of 
the driving light. The dashed-dotted line corresponds to the 
regularization by the two-body T-matrix, the dotted by the 
cut-off kc = Q.QlaQ , and the dashed line is the unregularized 
case. The solid line has A = 0. 

We studied the interaction of light with a two-species 
atomic superfiuid gas. The analysis of the quasi-particles 
followed the standard BCS theory [|lO|. We assumed a 
translationally invariant system. Atoms in a harmonic 
trap may be considered as locally homogeneous Q pro- 
vided that the trap length scale I = {Jijmijo)^!'^ is much 
larger than the correlation length, which for intraspecies 
\stvv ^ ^/kp and for the interspecies ^|| '^ ep/iAkp) 
|10| . Other notable assumptions were zero temperature 
and low atom density. In the future the present work 
could be extended to larger values of \ag\ and p by going 



beyond the BCS weak coupling limit and by including 
the cooperative optical linewidths and line shifts ||lj] . 
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